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Abstract: We find new black 3-brane solutions describing the "conifold gauge theory" 
at nonzero temperature and baryonic chemical potential. Of particular interest is the low- 
temperature limit where we find a new kind of weakly curved near-horizon geometry; it is a 
warped product A0IS2 xl 3 x T 1 ' 1 with warp factors that are powers of the logarithm of the 
AdS radius. Thus, our solution encodes a new type of emergent quantum near-criticality. 
We carry out some stability checks for our solutions. We also set up a consistent ansatz 
for baryonic black 2-branes of M-theory that are asymptotic to A0IS4 x Q ,1 > 1 . 
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1. Introduction 

Via the AdS/CFT correspondence Jj], ^[ |3|, electrically charged black holes in space-times 
with negative cosmological constant yield insights into the physics of strongly-interacting 
systems at nonzero density. Recall that the correspondence relates semi-classical gravity in 
d + 1 space-time dimensions to a strongly-interacting field theory in one fewer dimensions. 
There have been two kinds of approaches to these problems. In the bottom-up approach, 
a simple and phenomenological gravity model is constructed to which the AdS/CFT dic- 
tionary is then applied. It is assumed that there exists some field theory, of which we may 
have only a qualitative understanding, dual to this space-time, but using the correspon- 
dence, we can quickly and efficiently determine the phase structure, the equation of state, 
and transport coefficients. In the top-down approach, one considers a well-established 
AdS/CFT duality where the field theory is well-known but may have exotic symmetries 
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and field content. Given the complexity of the known dual pairs, calculations are often 
more difficult. However, they are worth the extra effort, because they allow us to make 
precise and reliable statements about actual strongly interacting field theories. 

In this paper, we take the second approach and construct a novel type of charged black 
hole (or, more precisely, black 3-brane) using the baryonic symmetry of the conifold gauge 
theory Q dual to the AdS§ x T 1,1 background of type IIB string theory. Our solution is 
similar to the Reissner-Nordstrom AdS§ black hole but is more complicated because the 
compact space T ' gets squashed by functions that depend on the radius. In the zero- 
temperature limit, we find that the near- horizon region becomes similar to AdS 2 xl 3 xT 1,1 
up to slowly varying logarithmic functions. The presence of the logarithms makes our IR 
solution a new kind of nearly conformal behavior. Thus, very interestingly, our solution 
exhibits "emergent quantum near-criticality," which could make it useful for exploring 
connections with condensed matter phenomena. 

Recently, there has been much interest in bottom-up approaches to study field theo- 
ries that undergo superfluid or superconducting phase transitions (see [||, |6| for reviews). 
Strong electron-electron interactions are believed to play an important role in the physics 
of high-temperature superconductors. Refs. M, ||, |j| |l0| use a bottom-up approach to 
model a strongly-interacting system of fermions at nonzero density, and find evidence for 
the existence of a Fermi surface. 

While the bottom-up approaches allow one to scan quickly through a number of simple 
gravity models and search for new phenomena, they have some disadvantages. A major 
issue is that the precise nature of the dual field theory is typically unclear, and one cannot 
be certain that it exists. Another disadvantage is related to the notion of a consistent 
truncation and its stability. The AdS/CFT correspondence in its original incarnation is a 
mapping between type IIB string theory in an AdS§ x S 5 background and M = 4 super 
Yang-Mills theory in 3+1 dimensions. In order to reduce a ten-dimensional string theory 
to a manageable five-dimensional gravity theory, a consistent truncation is made that elim- 
inates all but a small number of fields. The consistency of the truncation guarantees that a 
solution to the five-dimensional equations of motion for the remaining fields is also a solu- 
tion to the full ten-dimensional system. However, nothing guarantees that this solution is a 
global or local minimum of the action in the ten-dimensional setting; indeed, often it is not. 
The simple gravity models in these bottom-up approaches would be consistent truncations 
if fit into the larger AdS/CFT framework, and as such they may have instabilities. 

In top-down approaches, one of the easiest ways to charge a black hole is through the 
R-symmetry. The dual field theory often has an R-symmetry which maps to a gauge field in 
the gravity system. In the grand canonical ensemble, R-charge on the black hole translates 
into nonzero R-charge chemical potential [x in the field theory. In the most studied case of 
M = 4 super Yang-Mills in 3+1 dimensions, it is strongly suspected that at any nonzero 
/x, the theory is only metastable pi] , 12, 13|. Moreover, if \x becomes sufficiently large 
compared to the temperature, the theory becomes thermodynamically and perturbatively 
unstable as well [|J, ||, 

More generically, charged scalars, if their charge-to-mass ratio is sufficiently large, are 
a source of instabilities. Starting with refs. [17, 18], these charged scalars have been stud- 
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ied intensively in the context of modeling a superfluid or superconducting phase transition. 
Truncating such a charged scalar out of the gravity model also spuriously eliminates the 
phase transition. While it is true that the physical systems of interest often have a super- 
conducting phase transition at very low temperatures, from a theoretical point of view, it is 
of value to have a model where one may reliably go to low temperatures without worrying 
about such instabilities. 

In this paper, we construct a black 3-brane charged under a baryonic U(1)b symme- 
try that we hope is immune from the instabilities that plague R-charge black holes near 
extremality. We consider the well-known conifold gauge theory and its dual pair, type IIB 
string theory in an AdS 5 x T 1,1 background j§. The SU(N) x SU(N) field theory with 
bifundamental fields Ai and Bj, i,j = 1,2, has a global baryonic U(1)b symmetry. The 
corresponding U(1)b gauge field in AdS*, comes from the R-R 4- form with 3 indices along 
a topologically non-trivial 3-cycle. This realization of the U(l) symmetry makes our ap- 
proach different from the previous attempts to embed charged AdS black holes into string 
theory. In particular, the nature of the charged objects is quite different. The gauge invari- 
ant operators with baryonic charge in the conifold gauge theory have conformal dimensions 
of order N. The smallest such operator in the conifold theory involves determinants of the 
bifundamental matter fields. In the string dual, such an operator maps to a wrapped 
D3-brane which may be studied semi-classically. We are able to show explicitly that this 
wrapped D3-brane has a charge-to-mass ratio that is too small to produce an instability. 
This check, however, is insufficient to demonstrate the stability of our solution because 
one of the neutral fields may condense as the temperature is decreased. We demonstrate 
stability with respect to one seemingly dangerous neutral mode, but leave investigation of 
other modes for the future. We proceed with the optimistic assumption that our solutions 
are stable even at T = 0. 

In section ^, we review the details of the conifold gauge theory and its gravity dual. 
We also demonstrate a consistent truncation of type IIB supergravity (SUGRA) for the 
conifold background to a baryonic gauge field and two neutral scalars in five dimensions. 
Given the effective 5d Lagrangian, in section |3] we construct a metric, scalar, and gauge 
field ansatz for a baryonically charged black 3-brane that depends only on a single radial 
coordinate. The ansatz is invariant under a certain Z2 symmetry and leads to a system of 
non-linear ordinary differential equations. In section ||, we find their numerical solutions 
with AdS§ x T 1,1 boundary conditions at large r. Although we encounter difficulties at 
very low temperatures, the numerical work provides us with useful intuition concerning 
the T = solution. In section || we find the near-horizon series expansion in the T = 
limit and show that it has a nesx-AdSi structure. In section [6L we study the behavior 
of the smallest operator with baryonic charge in the conifold theory, the dibaryon. We 
show that its charge-to-mass ratio is too small to lead to an instability. In section ^ we 
carry out another stability check of our T = solution, this time with respect to a neutral 
Z2-odd perturbation. Because of nontrivial mixing with the U(1)r gauge field, the analysis 
requires us to generalize the ansatz of Section |2[ Section || contains some final remarks 
and discussion. In Appendix |A|, we find an analytic expression for our black hole in the 
small charge limit. In Appendix p] we consider black 2-branes carrying baryonic charges 
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in M-theory on AdS^ x Q ' ' , and develop a consistent ansatz and the equations satisfied 
by these solutions. 

2. Conifold Gauge Theory and Consistent Truncation 

The principal aim of this paper is to study the AdS/CFT duality in presence of a chemical 
potential for baryon number. In a general large N gauge theory, the operators carrying 
baryon number have dimensions of order N, which distinguishes them from the "mesonic" 
operators whose dimensions are of order 1 in the large N limit. In AdS/CFT, the objects 
dual to baryonic operators are D-branes or M-branes wrapped over non-trivial cycles of 
the internal manifold. In the maximally supersymmetric version, which relates the N = 4 
SYM theory to AdS 5 x S 5 , there are no baryonic operators. Their absence is related to the 
fact that S" 5 has no topologically non-trivial 3-cycles that could be wrapped by D3-branes. 
However, there are many known examples where the AdS/CFT duality relates AdS§ x Y, 
where Y is a Sasaki-Einstein manifold, to N = 1 superconformal gauge theories that have 
baryonic operators. The compact space Y typically has non-trivial 3-cycles, so that the 
wrapped D3-branes are topologically stable and the corresponding gauge theory possesses 
baryonic U(1)b symmetries. We would like to turn on a chemical potential for the baryon 
number in the gauge theory; in the string dual this translates into turning on the R-R 
4-form gauge field that couples to the wrapped D3-branes. 

We will present perhaps the simplest example of such a construction, in the context 
of the duality relating the AdS§ x T 1 ' 1 background of type IIB string theory to the M = 1 
superconformal "conifold gauge theory" 0. This duality is motivated by studying a stack 
of N D3-branes placed at the tip of the conifold, the Calabi-Yau cone zw — uv = 0. The 
explicit metric of its Sasaki-Einstein base, the T 1 ' 1 , is 

1 2 1 

ds?!,! = - ^2 (d9f + sin 2 6»j# 2 ) + - (dtp + cos M</>1 + cos 6> 2 # 2 ) 2 • (2.1) 

The conifold gauge theory with gauge group SU(N) x SU(N) is coupled to bi-fundamental 
chiral superfields Ai, i = 1,2 transforming as (N,N), and Bj, j = 1,2 transforming as 
(N, N). These superfields form doublets under the global symmetries SU (2) a and SU (2)b, 
respectively, and all of them carry R-charge 1/2. In addition to the U(1)r x SU(2)a x 
SU(2)b global symmetry, which on the gravity side is realized through isometries of T 1,1 , 
the gauge theory has a baryonic U(1)b symmetry under which 

A k -> e ie A k , B l e- ie B l . (2.2) 

The spectrum of gauge invariant operators splits into two sectors. The mesonic opera- 
tors, which are not charged under the U(1)b, have dimensions of order 1; the baryonic 
operators, which are charged under the U(1)b, have dimensions of order N. The low- 
est dimension examples of mesonic operators are Tr A{Bj of dimension 3/2, and TiAiAj, 
TrBiBj, Tr(AiAi — BjBj) of dimension 2. In general, the mesonic operators transform 
under the U(1)r x SU(2) x SU(2) geometric symmetry of T 1,1 , but are neutral under the 



- 4 - 



U(1)b- The lowest dimension operators carrying the U(1)b charge are, for example, det^4i 
or det^2- These are the m = ±N/2 states of the spin N/2 representation of SU(2)a- The 
general form of these dimension 3iV/4 operators, A m , m = —N, —N + l, . . . , N — 1, N, may 
be found in p9fl . They carry R-charge N/2, and we will normalize their baryon number 
to 1. The string theory objects dual to these operators are the D3-branes wrapping the 
(81, 0i, VO directions. Quantization of the (#2>02) collective coordinate gives rise to the 
N + l degenerate ground states corresponding to the chiral operators A m . 

Similarly, there exist chiral operators B m , m = —N, —N + 1, . . . , N — 1, N, such that 
Bn = det.Bi and B-n = det£>2- These operators have baryon number —1 and R-charge 
N/2; they are dual to D3-branes wrapping the (02,<p2,ip) directions. Replacing these D3- 
branes by anti-D3 branes we find objects of baryon number 1 and R-charge —N/2 that are 
dual to the antichiral operators B m that include detl?i- Since the U(1)r charge couples 
to mesonic operators and may lead to instabilities mentioned above, we will be interested 
in objects charged under the U(1)b symmetry only. The simplest such vertex operators 
are the (N + l) 2 products A mi B m2 which are dual to combinations of D3-branes wrapping 
both the (9i,(j>i,ijj) and (92,<p2,ip) directions. Turning on a chemical potential for U(1)b 
is expected to create a nonzero spatial density of such wrapped D3-branes. Our goal is to 
determine the background produced by them. We will use the simplifying assumption that 
the wrapped D3-branes are appropriately smeared over the T 1,1 coordinates, so that our 
solution will have the full SU (2) x SU(2) symmetry. 

The U(1)b gauge field in AdS§ x T 1,1 is contained in the components of the 4-form 



R-R-gauge field [ 2(| , C4 ~ A A W3 , where 



w 2 = g (sin 6*1^1 A d<pi — sin 62 d.82 A dfo) , 

^3 = 95 A uj 2 , ^ 2 - 3 ) 
55 = dip + cos 9id(pi + cos 6*2^02 . 

Our ansatz for the self-dual 5-form field strength will therefore be 1 

F 5 = -(F + *T) , 
9s 

2L 4 L 3 
T= —UJ2 Au 3 + —FAuj 3 , (2.4) 



4 _on„ /1 . L 



2 



Li 3\/2 

The normalization of the terms involving F has been chosen so that the kinetic term 
for F in the effective five-dimensional action is normalized canonically in the ultraviolet. 
Quantization of the 5-form flux requires that pit] 



L 4 = 4vr< 7s iV(a / ) 2 ^ . (2.5) 
lb 



x We define the Hodge dual of a p-form lo p in d space-time dimensions to be {*i^)i 1 i 2 ...i d _ p 
A. e J1J2 '" Jp ili2 id Uj 1 j 2 ...j p , where in a frame basis £12. ..d = 1. This is equivalent to requiring u> A *oj 
M 2 vol, where \w\ 2 = iff 1 "' 1 ■• • g^Ui l ...i p u j - s _... jv . 
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At non-linear order, we expect the additional components of F$ to cause a violation of 
the Poincare invariance of the 5d metric, and also to produce a squashing of the internal 
space T 1 ' . A minimal consistent truncation of type IIB supergravity that contains these 
effects turns out to be 



d s \ = e-l x ds 2 M + L 2 e x 



T) —4:7] 

j J2 [dOt + sin 2 6idtf) + e —gl 



(2.6) 



If x = f] = 0, this metric reduces to the direct product between the non-compact space 
M and T 1 ' . The scalar \ controls the overall size of T 1 ' 1 , while r\ introduces a stretching 
of the U(l) fiber relative to the two two-spheres. These scalar fields in AdS§ are dual to 
operators of conformal dimension 8 and 6, respectively. 

The above ansatz yields a consistent truncation of type IIB SUGRA with the effective 
five-dimensional lagrangian 

C eS = R-~e-lx+ 2 *F* v -5(d^) 2 -™(d flX ) 2 -V( V ,x) , (2-7) 
where the potential for the two neutral scalars is given by 

V(V, X) = ^e~f * + peri* (e" 6 " - Be"") . (2.8) 

The scalar kinetic terms and potential had been previously determined in [^]. However, 
the scalar coupling to the U(1)b gauge fields was not considered there. Indeed, fl2.7| ) shows 
that the gauge kinetic term depends on the scalars; to study the baryonic black holes we 
need to include the squashing of T 1,1 . 



A seven-dimensional analogue of T ' is the Sasaki-Einstein space Q ' ' [23]. The 
background AdS± x Q ' ' has two kinds of baryonic symmetries corresponding to two 
topologically non-trivial 5-cycles that can be wrapped by M5-branes. In Appendix [b] 
we consider black 2-branes carrying one of these baryonic charges, and find a consistent 
truncation which includes a squashing of Q ' > . 



3. Equations of Motion 

Using ( |2.7D , we will look for time-independent charged black 3-brane solutions. We impose 
rotation and translation symmetry in the 3 spatial directions, but the Poincare symmetry 
is obviously broken. Therefore we will use the ansatz 

ds 2 M = -ge- w dt 2 + — + J2 5>^) 2 > ( 3 - X ) 

where g and w are functions of r. This choice of parametrization (see, for example, 
will prove useful in simplifying the form of the equations of motion. To turn on the baryonic 
charge density and chemical potential, we need to consider only the time component of the 
U(1)b gauge field, A = <S>(r)dt, so that the field strength is F = dA = <S>'dr A dt. We will 
also assume that the scalars x an d f] depend on the radial coordinate r only. 
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We should note that our full ten-dimensional ansatz fl2,4| )-(l2,q) preserves the Z2 space- 
time inversion symmetry where (t, x) —¥ (— t, —x) is accompanied by the interchange of the 
two 2-spheres, (#1, 4>\) -H- (62, 4>2)- The forms 0J3 and 0J2 change sign under this transforma- 
tion, but the terms dt A dr A 0J3 and dx 1 A dx 2 A (fx 3 A 0J2 present in F§ are invariant. In the 
gauge theory, this Z2 symmetry appears to correspond to (t, x) — > (— t, —x) accompanied 
by Ai o Bi. 

With this ansatz, the effective one-dimensional lagrangian is 

L cS = -5rV~H?/ 2 + \x 2 ) + VeK 2 ""^' 2 - 3( 5 r 2 )' e -^ - 
The equations of motion following from this lagrangian are 

,11 , ,jf s , 5' , 5r^ /2 , lOr , 2 ^ *%+2„-§ x __3_^ 



r g 3 9 / Wg 20c? % 

./^ 3 , 5' , 5r^ /2 , lOr /2 \, , &\ w+2v _i x _ J_dV 



r 3 e~^ w V . 


(3.2) 


= , 


(3.3a) 


= , 


(3.3b) 


= , 


(3.3c) 


= , 


(3.3d) 


= , 


(3.3e) 



rf' + rf [- + - + — rf* + —X*\+—e 

r g 3 9 / 10g lOg ary 

\r 2 3 
W ' + ^(3r/ 2 + 2 X /2 ) 

where primes denote derivatives with respect to r, and V is the scalar potential defined in 

& 

Equation ( |3.3d| ) can be integrated once yielding 

$' = ^ e -^-2i7+|x j (3.3f) 

where Q is an integration constant related to the charge of the black hole. This is the 
conservation equation for baryonic charge. From now on, we will use ( [3.3f] ) instead of 
( 3.3d ). Plugging ( |3.3f| ) into equations ( p.3a| ) ( |373c ), we get the added bonus of eliminating 



the dependence on w, so the remaining equations are 



.// , . ./ ( 3 , 9' , 5r /2 lOr /2 



V" + V' - + - + ^ + ^ + 77f7T^ +f « - = , (3.4b) 

9' + 9 (~ + f r/ 2 + + ^ e- W + 3^ = 0- ( 3 ^) 

To reduce the system from five to three coupled differential equations was the main moti- 
vation for using the ansatz (^ 



Q 2 




3 ay 


Wr 6 g 




20g d X 


Q 2 




1 ay 


Wr 6 g 


10g dr] 


/2 ) + 


6r 5 


3 



x" + x' - + - + ^v' 2 + ^x' 2 - -rhr^ 2 ^ x -~ = o, (3.4a) 

\r g 3 9 
// , / / 3 , 9' , 5r /2 lOr , 2 
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4. Numerical Solutions 



Solving the coupled non-linear equations ( |3.3a )-( 3.3f ) is in general a difficult task. We 



will mostly rely on numerical work, but in some limits will be able to present analytical 
formulae. The simplest situation is when Q = and we find the well-known black 3-brane 
solution in AdS$. In this solution, the shape of T 1 ' 1 does not depend on r, i.e. the scalars x 
and r] vanish. For small values of Q we can use perturbation theory in this small parameter 
to obtain an analytic expansion of the solution. This exercise is carried out in Appendix [A| 
where we find that the scalars x an d f] are now of order Q 2 and acquire a dependence on r. 
In the next section, we present some analytical results in the extremal limit. However, for 
intermediate values of Q, we know of no good analytical methods and resort to numerical 
ones. 

4.1 Setup for Numerics 

Finite-temperature solutions are found numerically by a standard shooting technique. The 
numerical solver is seeded close to the boundary, which is located at r — > oo, by using a 
series expansion that imposes the correct boundary conditions. 

The first task is to determine the boundary conditions. All our ten-dimensional metrics 
must asymptote to AdS§ x T ' at large r. This means that the asymptotic boundary 
conditions that we require as r — > oo are 

2 

^0, g = r - + 0(L 2 /r 2 ) , (4i) 
ri , x -> . 

In order to describe states in the dual field theory at nonzero baryon chemical potential, we 
also require $ — > <3?o> f° r a constant <l?o that will be related to the chemical potential shortly. 
Generically, solutions satisfying these boundary conditions will have an event horizon at 
some value r = where g(r) vanishes. The standard boundary conditions required at the 
horizon are that w, 77, and x should be finite. 

From ( 3,3a| )-( |33f} ) one can work out a series solution at large r that satisfies the 



boundary conditions ( |4.1[ ): 

w = 0((L/r) 12 log 2 (r/L)) , 

5 = 5 + ^ ! + T^ + O((W 10 log 2 (r/^)) , 

^ = ^-^ + 0{{L/rf\og(r/L)) , (42) 



200r° r° 



ri 10 ) 



^ o^y ) + ^ +0((L/r)10log(r/L)) 

All higher order terms in the series are determined in terms of g2, X8> an d Q- 

To proceed further, it is useful to review some of the symmetries of our ansatz. The 
equations of motion ( |3.3|) and the boundary conditions are invariant under some scaling 
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symmetries which act with the charges summarized in table |]. We say that a quantity X 
has charge q under a scaling symmetry if 



X -> \ q X 



(4.3) 



The first symmetry in table |T] is a formal way of expressing the arbitrariness of a choice 



Symmetry 


e~ w 


9 


$ 


V 


X 


t 


X 


r 


L 


type A 

















1 


1 


1 


1 


type B 





2 


1 








-1 


-1 


1 






Table 1: Charges under the scaling symmetries of the equations of motion ( B.3) and of the boundary 



conditions (4.1). These charges are defined as in (4.3). 



of units in the bulk. It can be used to set L = 1. The second symmetry can be used to 
put r/j = 1, but when shooting from the boundary it is more useful to employ the same 
symmetry to set g<i = —1 instead. We can also set $o = 0, by using the fact that the 
potential $ appears only through its derivatives in the equations of motion, and does not 
appear in the metric. Having eliminated <£o an d fixed Q, we are left with two parameters, 
Xs and tjq, that need to be tuned in order to match to regular solutions at the horizon. 
Trial values for these parameters are determined from the small Q (large T) expansion of 
Appendix A, and the results are used to go to progressively smaller temperatures. 

Thermodynamic quantities such as the energy density e, entropy density s, temperature 
T, charge density p, and chemical potential /i, can be computed from the following formulae: 



3 92 



2k\L 
Q 



2k 2 5 L 2 



ft 



2vrr| 

$0 - $h 
L 



g'(r h )e 



-w h /2 



4-7T 



(4.4) 



where K5 is the five-dimensional gravitational constant, = io(r/j), and $/j = ^(r^). 
Using the parameters from [21] we find 



1 



27iV 2 



k\ 64vt 2 L 3 ' 
The energy density can also be computed from 

e = 4 ( Ts + ^ 



(4.5) 



(4.6) 



which follows from e = Ts — p + pp and the tracelessness of the stress tensor, e = 3p. 



4.2 Results 

Using the shooting method described above, we were able to find numerical black hole 
solutions for a fairly large range of temperatures. Here are the main features we observe: 
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Figure 1: The behavior of the horizon values of the scalars as a function of T j \i. The dotted line 
is the high-temperature behavior obtained from the low Q expansion in appendix ^[ According to 
the zero-temperature expansion developed in section ^, r/h and \h should diverge at T = 0. We 
believe we are far from this regime. 



Our numerical results agree with the analytical computations presented in Appendix |A] 
in the limit of high temperatures. See figure |l] for a comparison of the values of the 
scalars at the horizon found numerically to those predicted by the analytic formula 



• We were able to construct numerical black hole solutions only for temperatures higher 
than T ~ 0.0005/U because of loss of numerical precision at lower temperatures. We 
believe the lowest temperatures we attained are not low enough to provide a thorough 
check of the analytical zero-temperature horizon expansion constructed in the next 
section. In fact, we will argue towards the end of the next section that we expect the 
zero-temperature expansion to become a good approximation to the near-extremal 
solutions when loglog/x/T 3> 1, which is beyond the range of temperatures where 
our numerics are reliable. 

• As one decreases the temperature, the Bekenstein-Hawking entropy seems to ap- 
proach a non-zero value (see figure ||) . This fact will be confirmed analytically in the 
next section, where we find that as T approaches zero, s/p ~ 2.09. However, as will be 
discussed in the next section, there are stringy effects that make our solution trust- 
worthy only down to an exponentially low temperature of order /Ue - const - x (s , sA0 2,/3 - 



• Lastly, one may worry that at low enough temperatures the curvature of these back- 
grounds might get large and the supergravity approximation might break down. As 
can be seen from figure ||, the horizon value of the ten-dimensional Riemann tensor 
squared is uniformly bounded from above. In fact, in the following section we will 
show that all curvature invariants should vanish at the extremal horizon. 

5. Near-AdS 2 Near Horizon 

We now find a horizon series expansion at zero temperature. We will restrict ourselves to 
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Figure 2: Entropy density normalized by charge density, as a function of temperature over chemical 
potential. The entropy is seen to go to a nonzero constant as T — ¥ 0, in agreement with the zero- 
temperature expansion discussed in section ^. 
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Figure 3: The five-dimensional Ricci scalar (left), and the square of the ten-dimensional curvature 
(right) at the horizon, in units where L = 1, as a function of T//z. As discussed in section |[ all 
curvature invariants evaluated at the extremal horizon vanish. 



the set of equations fl3.4|), since one can always use ( |3.3ej ) and ( |3.3f ) to find w(r) and <3?(r) 
afterwards. In the following, we use the symmetries in Table |l| to set L = = 1. 

A guess for the zero-temperature value of the charge Q that appears explicitly in 



equations (3.4) can be found from the following line of reasoning based on properties of 
the nonzero temperature solutions. At the horizon, g(fh = 1) = 0; since g must be positive 



outside the horizon, we need g'(l) > 0. Evaluating eq. ( 3.3c ) at the horizon, we get 
1 



</(!) 
which implies 



6 



(I6e 6% + 8e 4 * h - 48e 5llh+iXh + Q 2 e 4rih+8xh ) > , (5.1) 



bi < e 4xh < b 2 , 



where 



n,2 



■1 + 6e 5rih ± y/1 - 12e 5 ^ - (Q 2 - 36)e 10 ^) 



(5.2) 



(5.3) 
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For Q > 6, both b± and 62 are negative or complex for any 77^, resulting in no possible 
range of Xh- F° r Q < 6, the smallest positive value of rjh for which b\ t 2 are rea l is the one 
for which b\ = 62. The fact that there are no solutions for Q > 6 implies from ( |4.4j ) that 



s 4ir 2i , 
_ = _>_ K2 .09, (5.4) 

an inequality that should hold at all T within the supergravity approximation. 
When b\ = 62 5 we have 

r h = ~]og(6-Q) , Xc = --Uog(6-Q)-ilog^ , (5.5) 
5 20 4 4 

which implies that r/ c and Xc go to infinity as Q — > 6, while 

Q 3 

t/ c - 4x c = log — -> log - as Q -»• 6. (5.6) 

A reasonable guess is that at zero temperature Q = 6 and the scalars x an d f] diverge at 
the extremal horizon, while r\ — 4^ approaches log |. 

With Q = 6, we use the method of dominant balance to find the asymptotic behavior 
of zero-temperature solutions at the horizon. By adding arbitrary power series to these 
dominant terms, we find we can satisfy the equations of motion. We obtain a solution of 
the form 

1 , /2187 \ 1063 

y = log r H r + • • • , 

A 20 & V 16 / 1000 

1, , . 463 

V = --log(18r) + — r -\ , 

5 250 , g ^ 

/ 93312 • (12) 1 / 3 



9 = f 13 / 3 



V 25 



5 77, _ „ 

w = 1 log r + w + 

36f 30 & 



where 



1 , (5.8) 



and the dots stand for regular Taylor series. 

As f —> 0, goo = —g e~ w has an essential singularity. This prompts us to introduce a 
better coordinate, 

y = 3 -193/120 2 -107/120 5 3/40 e «;/2 g 5/(72f) ^ 

which becomes large near the horizon. In terms of this coordinate, and reinstating the 
factors of L 2 , the near- horizon metric becomes 

9 ( /1 \ 37 /O.D 1 9 1 /-i \ 1 I A. 1 o\ T O. ( 1215 

as 



'10 — ~T 

yZ 



+ 



(logy)- 37 / 20 ^ 2 + ^yi7l(logy) 1/4 V) + (^) ' (logy)-^dx'- 

) (log y) 1/4 £ + sin2 W?) + Y ( If ) ( lo s »)" 

/ A — 1 \ / 



fh 
(5.10) 
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This is a warped product AdS2 x M 3 x T 1,1 with warp factors that are powers of the 
logarithm of the AdS radius y. The appearance of the logarithmic warp factors makes it a 
novel type of "nearly conformal" IR behavior. Note that the 3d spatial components of the 
metric tend to zero in the IR as g^j ~ (log y)~ 1//12 <%; this is an important difference from the 
RNAdS metric where they approach a constant. Remarkably, the extra logarithms actually 
reduce the curvature, so that we can trust the supergravity approximation everywhere: 
we have checked that the lOd Kretschmann invariant and the 5d Ricci scalar scale as 
RabcdR abcd ~ (logy) -1 / 2 and R ~ (logy) -1 / 3 respectively. Such an asymptotic reduction 
of curvature due to the appearance of logarithms also occurs in the UV region of the warped 



deformed conifold [22, 25]. 

One can use the asymptotic solution to estimate what happens at temperatures so low 
that the horizon is located deep inside the near-AdS^ region, at <C 1. If we assume 
that ( [5.7D approximately holds down to the horizon, we can estimate using fl4.4| ) that 

_ 5 

T//x ~ e 72f h up to power law corrections in fh and numerical factors. From ( |5.7[) , one 
can then compute the values of the scalars at the near-extremal horizon: 

rih = *Xh = — log tv, + O(r h ) = ~ log log £ + 0((T/^) ) . (5.11) 

We therefore expect the approach to the extremal limit to be very slow and hard to inves- 
tigate numerically because one has to go to exponentially small temperatures. 

This numerical suggestion of the existence of an exponentially small scale compared to 
the chemical potential [i in the conifold gauge theory is corroborated by a stringy argument. 
In the type IIB string theory context, our solution cannot be trusted for arbitrarily low 
temperatures. 2 In the T = solution the ^-circle shrinks to zero size at the horizon (see 
the last term of (|5 . 1 0|) ) , and the standard approach is to T-dualize along this direction for 
logy > const, x (g s N) 2 / 3 where the size of the circle becomes of the order of the string 
length y/c/. Estimating the temperature at which the size of the ■(/'-circle at the horizon 
reaches the string scale, we find 

log|~( 5s iV) 2 / 3 . (5.12) 

It is remarkable that application of string theoretic arguments to our black brane 
suggests that in the conifold gauge theory at baryon chemical potential /x there exists 
an energy scale of order ^ e - const - x (9sN) 2 3 F rom the point of view of the type IIB a- 
model with coupling ~ (g s A)~ 1//4 , such a scale can arise only non-perturbatively. Below 
this exponentially small scale, the gauge theory presumably exhibits some new effects 
that can be studied by T-dualizing the solution and lifting it to M-theory. While the 
background ( |5.7| ) has a non-vanishing Bekenstein-Hawking entropy and is smooth, stringy 
effects become important when the -0-circle becomes small, and these effects could remove 
the conflict with the conventional statement of the Third Law of Thermodynamics | )26| ]. 
Since the type IIB backgrounds can be trusted when the ^-circle is above the string scale, 
we expect the entropy to be approximately constant for ^ ie ~ const x (9s Ar ) 2/ ' 3 < T <C /i, i.e. 



2 We are very grateful to J. Maldacena and E. Witten for pointing this out to us. 
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when the horizon is inside the near- AciS^ throat. It would be very interesting to provide a 
microscopic origin of this entropy of order N 2 by studying the underlying D3-brane system, 
which involves two intersecting stacks of D3-branes wrapped over the T > . We hope to 
return to these issues in the future. 

In order to match the asymptotic solution ( f>.7\) to the appropriate behavior at the 
boundary, we generically need a two-parameter family of solutions, in addition to the 
integration constant w, which is related to a choice of units for time in the boundary 



theory. Following [24], we linearize around this leading behavior, to find non-analytic 
pieces that we might have missed. There are two solutions to the linearized equations, of 
the form 

8 x = e <*/?r ci {l + ---) , 

S v = 45 X + e a/f f C2 (l + ■■■), (5.13) 
Sg = e a / f r C3 (l + ---) , 

where the corrections are just series in positive powers of f. The values of a and Cj for the 
two solutions are 

" = ' Cl = 20 ' c 2 = ci + 1 , c 3 = ci + — , (5.14) 

and 

5 , r— ^ 179V2T- 189 13 



a = ~TuS^ > ' Cl = 360 ' C2 = Cl ' c 3 = ci + y- \P-H>) 

For a scalar in AdS%, the behavior of linearized modes as a function of the radius y is 

cf>(y, t) ~ 0o(t)y 1_AlR + A(t)y A ™ , (5.16) 

where </>o(£) is the source for an operator of dimension Air, and A(t) is its expectation 
value |2(|. In the case at hand, since the background is AdSz up to powers of logy, we 
expect ( 5.16j ) to hold up to similar powers. The perturbations ( 5. 13] ) are indeed of this 



form. The first of these solutions corresponds to a source for an operator of Air = 2; the 
second to a source for Air = (1 + v21)/2. Both of these operators are irrelevant in the 
nearly conformal quantum mechanics found in the IR; hence, they produce perturbations 
that decay at large y. 

Another type of perturbation that we can study easily is a minimally coupled massless 
scalar, for example the dilaton. In the near-^ldS^ region, the minimal scalar solutions 

behave as a + b y (log y) 4//5 -I . We identify the dual operator as a marginal operator with 

Air = 1. In the case of the dilaton, the operator is exactly marginal because the string 
coupling is a parameter in our solution. 

The extremal solution described above is reminiscent of the "run-away" attractor flows 



described in [27]. Typically, an attractor flow is a set of solutions to the supergravity 
equations of motion where the ultraviolet behavior is not universal, but in the infrared the 
scalars approach fixed values. In the run-away case at least one of these fixed values is 
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at infinity. Our extremal solution would be of run-away type because the scalars 77 and 
X diverge at the extremal horizon, even though the combination 77 — 4\ stays finite. An 
important difference between the extremal solutions from 27] and the one we found is that 
in [27] the entropy density vanishes at extremality while in our case it does not. 



6. Baryonic Operators and D3-brane Probes 

In this section we investigate whether the baryonic black branes constructed above are 
stable with respect to condensation of baryonic operators in the conifold gauge theory. 
As was already described, these operators are dual to wrapped D3-branes in the dual 
supergravity. These wrapped D3-branes act like charged particles; there is a competition 
between the gravitational attraction and electrostatic repulsion between the particle and 
the charged black brane. To test stability, we propose a simple thought experiment. If we 
place such a wrapped D3-brane into the geometry and the brane falls into the black hole, we 
conclude the black hole is stable. However, if the brane finds some meta-stable minimum 
outside the horizon, we conclude the black hole is unstable and more such wrapped D3- 
branes can bubble off the horizon, find their way to the minimum in the potential, and 
reduce the baryonic charge on the black hole. It is possible that more exotic bound states 
of D-branes may lead to instabilities. We leave an investigation of such issues to future 
work. 

6.1 Charged Particles in Reissner-Nordstrom AdS 

Instead of beginning with the D3-brane case, let us first consider the simpler case of a 
particle of mass m and charge q. The action for such a particle is 



m I ds 



dx» 



ds 



q I dsAf.-j- 



dx» 
ds 



We restrict to a metric of the form ( ]3.1| ) with a radial electric potential <E>(r) 
a particle sitting at some fixed r and x, the action becomes 



dt 



me 



-w(r)/2 



V§F) + q®{r) 



(6.1) 
A t {r). For 

(6.2) 



in the gauge s = t. Equation ( |6.2[ ) shows that the potential for this particle as a function 
of r is given by 



V{r) = me"" ,(r)/2 v / 9W + g*(r) . 



(6.3) 



The first term is the gravitational attraction of the black hole, and the second term the 
electrostatic repulsion. 

To proceed further, we need explicit results for w, g, and <£. We will return to the 
baryonic black hole momentarily, but let us first consider what happens for the simpler 
case of a Reissner-Nordstrom black hole with negative cosmological constant. Given the 
Einstein-Maxwell action 



Sem 



1 

2^2 



■d+l. 



-9 



R + 



d(d - 1) 
L 2 



1 

4? 



d d+1 Xy 



(6.4) 
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there is a black-bran e solution with a horizon at r = and with w = 0, 

» = S( 1 + ( 3 2 (7) 2 "-( I + < ? 2 >(7) d ) ' < 6 ' 5 > 



. i rh d — 1 



k L V d — 2 
The corresponding Hawking temperature is 



r / 



1 



(6.6) 



= 4^ • (6 - 7) 

Without loss of generality, let's assume Q > 0. The RNAdS black hole becomes extremal 
at Q = Q max = yj where T# vanishes. 

Using these formulae for g and <£, one can check that for values of q/m < q^it/m, V(r) 
is increasing monotonically for all values of Q smaller than <3max> while for large values 
q/m > q C rit/ m i V(r) has a minimum at some r = r*(Q) provided that the charge Q of the 
black hole is larger than some critical value Q c that depends on q/m. An expression for 
the critical charge <7 cr ; t can be found by requiring that the minimum of V(r) occurs right 
at r = rh when Q = Q max - The equation V'(rh) = is easily solvable in this case, yielding 

fet = m- , (6.8) 

e 

independent of the number of dimensions of the RNAdS space. 3 

This picture of an instability is essentially the classical limit of the superfluid or su- 
perconducting instability studied in refs. 0, In these papers, the charged particle is 
replaced with a charged scalar field ip. One solution to the equations of motion is tp = 0, 
but for Q > Qc, there exists a second, more stable branch of solutions with tp 0. The 
classical analog of this second branch is a cloud of charged particles sitting at the minimum 
in V(r) described above. Moreover, for large m and q, the classical and field theory results 
for Q c agree. 

6.2 The Wrapped D3-brane 

The probe D3-brane action takes the form 

Sds = -M3 J d^e-V-detGab + e/xg J C 4 , (6.9) 

where ^ = (27r)~ 3 (a / ) -2 , a' is related to the string tension, and <f> is the dilaton. For our 
supergravity solution, the dilaton is a constant we take to be related to the string coupling, 
g s = . The parameter e is equal to one, but we leave it arbitrary so that we can tune the 
charge of the D3-brane. We assume the metric ansatz ( |2.6j ) and fl3.1|). 



3 This ratio has some broader significance, as emphasized by [Eq] in the context of a weak gravity bound. 
If extremal Reissner-Nordstrom black holes are to be able to decay, there must always exist at least one 
particle whose charge-to-mass ratio is greater than n/e. 
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We assume the probe brane sits at constant 82, (fa, %i, and r and wraps the remaining 
four directions, including time. From £03 and the ansatz for F§ fl2.4|) we deduce a potential 
for the D3-brane: 



V{r) 



3N 



e— — — $(r) . 



(6.10) 



4L v ^ w 8L 

Comparing (6.10) with ( |6.3[ ) , we see that in the UV where the scalars r/ and x are 
negligible, e w ^ ~ 1 and 5 ~ r 2 /L 2 , the wrapped D3-brane corresponds to a particle in 
RNAdS with m/q = y/2. As shown at the end of section |6.1| , this ratio corresponds to the 
critical value m/q cr it at which in RNAdS the potential V(r) has a local minimum at r = 
in the extremal case. There is no reason to expect that the critical value of m/t/ cr it should 
be the same in the presence of scalars, but we will now show that this is nevertheless true. 

The simplest demonstration is to plot V(r) numerically. At the lowest temperatures 
we can access, there is no minimum in the potential when e = 1. However, for e > 1 and T 
sufficiently low, there is such a minimum, suggesting the D3-brane does indeed have this 
critical ratio of charge to mass. Indeed, the minimum is observed to occur for e > eo, where 
eo is some critical value larger than 1. In figure ||, we plotted the dependence of this critical 
value on T/fi. We see that as the temperature goes to zero, the critical value goes to 1. 
The position of the minimum moves towards the horizon, showing that wrapped D3-branes 
have the marginal ratio of charge to mass that means they barely escape condensation. 




1.0 



T//i 



Figure 4: Critical value of the parameter e above which condensation occurs, as a function of T/ \x. 
Just above eo, the position of the potential minimum goes towards the horizon as T — > 0. 

Some analytic support for this claim comes from studying the conditions that the 
minimum in the potential disappear, V'(r*) = V"(r#) = 0. Restoring e = 1, these two 
vanishing conditions imply the relation 



(6.11) 



However, we know this condition is satisfied in the T — > limit at the horizon from 
fl5.6| ). Thus the system just begins to become unstable at T = 0, and we expect no phase 
transition. 
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More explicitly, we examine the potential ( 6.1C| ) in the T — > limit. From the zero- 



temperature series expansion (p77|), we find that 

5 i - / 394 • 1 S 1 / 6 \ 
V(r) = iVe _ W^~^( r _ i)i/20 (i- e )_ ^^(r - 1) + 0(r - l) 2 , (6.12) 

where we set L = = 1. Thus, the leading term in this series expansion vanishes when 
e = 1. 

7. Another Stability Check 

The results of the previous section show that the black holes with baryonic charge we 
have constructed are stable against the simplest condensing operators that carry nonzero 
baryonic charge. Recall that all such operators have large dimensions of order N . One 
may worry about a different kind of instability where at low enough temperatures there is 
a phase transition driven by operators that are uncharged under the baryonic symmetry. 



That such a phase transition is in principle possible was noted in [29, |30|, |31J for the case of 
the RNAdS black hole. In that case, all uncharged operators with UV conformal dimension 
smaller than 3 (when the gauge theory is 3+1-dimensional) could trigger such an instability. 

In this section we will study the stability of certain modes in the IR neai-AdS2 region 
of the T = solution. For a minimally coupled scalar, 



^± = \±\\ + ^L\ dS2 . (7.1) 



2 V 4 

Hence, the Breitenlohner-Freedman (BF) stability bound in AdS2 is 



m 2 L 2 > -- , (7.2) 



1 

J AdS 2 ^ 

and the dimension Afr_ is allowed only in a narrow range of m 2 above this bound p0| . 
When the BF bound is violated the dimensions become complex, and the modes exhibit the 
oscillatory behavior as a function of the radius that is characteristic of "bad tachyons" (see 



for example [32]). This simple analysis is not directly applicable to the case of interest to 
us because the background is not exactly AdS2, and the scalars are not minimally coupled. 
However, we will adopt a similar stability criterion. In the linearized approximation, an 
instability will be associated with the presence of complex dimensions and the ensuing 
oscillatory behavior of modes. 

The baryonic black 3-branes studied in this paper come from a ten-dimensional type IIB 
construction, so one can study fluctuations of various (uncharged) supergravity fields and 
check their stability. A full analysis of all possible modes is beyond the scope of this paper. 
However, we will demonstrate stability with respect to perturbations associated with a field 
theory operator of UV conformal dimension A = 2. This operator is Tr(Ajylj — BjBj), 
and its dual supergravity field is the "resolution mode" of the conifold, A, that allows the 
two S 2l s to have different sizes. This is the most relevant mode that is odd under the Z2 
space-time inversion symmetry accompanied by the interchange of the two 2-spheres. We 
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thought that this mode was the most likely to cause an instability because it saturates the 
BF stability bound in AdS 5 x T 1 ' 1 ; luckily, as we show, it does not destroy the stability of 
the near-Ac^ solution. 

The challenge here is that A mixes with the time components of certain gauge fields, 
already at the linearized level. We were nevertheless able to find a consistent set of super- 
gravity fields that include A and that decouple from all other fluctuations, providing a more 
general (non-linear) consistent truncation of type IIB supergravity than the one considered 
in section |2|. Indeed, the consistent truncation we find reduces to the one in section |in a 
particular limit. The full ten-dimensional ansatz and the effective five-dimensional action 



are given in section |71|. In order to examine the stability of the near-AciS^ geometry, in sec- 
tion [7^ we linearize the equations of motion around the baryonic black brane background 
and develop a horizon series expansion at zero temperature using the explicit solution dis- 
cussed in section [5[ Remarkably, the linearized equations contain mixings between the 
modes and effective mass terms that stabilize all potentially unstable modes within this 
ansatz. 



7.1 A More General Consistent Truncation 

A consistent truncation that extends (|2.4|)-( |2~6| ) to include the resolution mode of T ' can 
be constructed as follows. Compared to ( |2.4[) - (|2.6D , this truncation has three additional 
fields: the scalar field A, which is the resolution mode of the conifold, and the spin-one 
fields Ar and Ar, which mix to give a gauge field corresponding to the R-symmetry of the 
gauge theory as well as a massive spin-one field. The metric ansatz is 

a ri+\ 



Defining 



dsj 



e~ 5x/3 ds 2 M + L 2 e x 



6 



-(d0f + sm 2 #i# 2 ) 



+ 



-(del + sin 2 9 2 



-iii 



+ 



95 + 



V2L 



A 



R 



2 i 



(7.3) 



55 = 95 + 



V2L 



Ar , 



the self-dual five-form can be written as 
1 



^5 

T 



(F + *F) 

9s 
2L 4 

— uj 2 Au 2 Ag 5 + ■ 7rj =F Aw 2 Ag 5 A 
L 3 



9\/2 J 



+ 



—e s *vo1m H i=<i 3 

L 3\/2 

L 2 

-.e s 



Ar A dg 5 A dg 5 , 
cosh(2A) *ju F — sinh(2A) *m Fr 



(7.4) 



(7.5) 



A UJ2 



+ 



6V2 
2V2 



cosh(2A) *m Fr — sinh(2A) *m F 



Adg 5 



-4x-4»7 



* M (Ar + A R ) A g£ , 
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where we have defined the field strengths F = dA, F R = dA R , and F R = dA R . 

The effective five-dimensional action for this consistent truncation can be written as a 
sum of a bulk piece and a Chern-Simons term: 



S 

The bulk lagrangian C is given by 
10 



J d 5 x^C + S CS ■ (7.6) 



C = R- y (5 mX ) 2 - Kd^V? ~ (d^f ~ V(rj, X , A) 



cosh(2A) (F^F^ + FfrFg) - 2 sinh(2A)F^F^ 
4 

J? 



4 

\*****F$E% - ^e- 4 ^« + A*f 



(7.7) 



where 



V(t], X ,X) = -^e-f x + -^e-f x (e- 6r 'cosh(2A)-6e- r 'coshA) . (7.8) 
The Chern-Simons part of the action is 

Scs = ^Ja r AF r AF r -^=JAAFAF r . (7.9) 

From now on, we restrict to a time-independent background where we have rotation 
and translation symmetry in the three non-compact directions x l . The most general ansatz 
with these symmetries is where the scalars x> V> an d A depend only on r, and where 

A = <Z>{r)dt , A R = <S> R (r)dt , A R = $ R (r)dt , 

generalizing the ansatz used in section ||| In this case, the equations of motion following 
from the effective action (|7.6| ) admit two conserved charges Qb and Q R associated to 
the baryonic symmetry and to the R-symmetry of the gauge theory, respectively. The 
conservation equations take the form 



$'cosh(2A) - l>' /? sinh(2A) = 9^. e -^-^x ? ( 7 lla ) 



$' R cosh(2A) - $'sinh(2A) - ^ e ~ Wx & R = 9*. e -fr»-*t+fr ; ( 711b ) 

where, as usual, primes denote derivatives with respect to r. The first of these two equations 
is a generalization of ( 3.31 ) . In the UV, the scalars are negligible and the above two 
equations reduce to <fr' ^ ^Qb and & R — \<&' R ~ ^Qr, justifying the interpretation of Qb 
as the baryonic charge and of Q R as the R-charge. 
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7.2 Horizon Expansion of Linearized Fluctuations 

As mentioned in section ||, our background is invariant under a Z2 symmetry that acts 
by flipping the sign of the non-compact coordinates, (t,x) — > (— t, —x), and interchanging 
the two spheres in the compact space, (9\,(f)i) f-> (62, 02 )• Recall that in the background 
solution only x, rj, and <1> are nonzero. The fluctuations around this background are 
distinguished by their parity properties under the Z2 symmetry: 5&, 5x, and 5rj are even, 
while 5& R , 5<t R , and 5X are odd. Due to the Z2 symmetry of the background, the even 
and odd linearized fluctuations cannot mix. Here we are interested in the resolution mode 
5 A, so we will focus on the mixing among the odd fluctuations. 
The linearized equations are 



3 \ A 

e2 w J L J L 

5 _( r 3 e -4r,+ |x+^^)' _ l| e -^-4x(5$ fl + S $ R ) = , (7.12) 



^( r 3 e 2r?-|x+^^)' _ 2ge 2r i-U5\'& - -^e" 4r '- 4x (5$ iJ + 6$ R ) = 



r 3 e X K ' * L 2 



where x, 77, 5, and ui are evaluated at their background values given in ( |3.31| ) and 
for the zero-temperature extremal solution. Let's focus on this extremal solution and find 
a series expansion in r — 1 (we set L = = 1). This calculation is similar to that of the 
non- analytic contributions to the background given at the end of section ||. Since ( 7,12; ) 



is a system of three second order differential equations, there are six linearly independent 
solutions whose leading behaviors are of the form 

5$ R = e Tr bl ( 5<Z> { y + StoQr + ■■■ 



, R = err \ov R -t- o\e R 



S* R = eTf>» (6*%> + 8*%r + ...) , (7.13) 
5X = e ^eri?f 6 3™i ( 6X (o) + JA (i) f + . . 



with f = r — 1 as in (|5^). The coefficients a and 6j are given in table ^. All six solutions 
satisfy the U(1) R charge conservation condition ( [7.11b ) at the linearized level. Of the 



modes that do not grow with y (I, III, IV, and V), only mode III requires a non- vanishing 
Q R , while the others satisfy the conservation equation with Q R = 0. The crucial fact is 
that a is real for all six solutions, so there are no oscillatory solutions in the near-Ad^ 
region. The absence of oscillatory solutions means that the black 3-branes with baryonic 
charge constructed in the previous sections are likely to be stable with respect to the 



perturbations ( 7.12 ). 

We discussed in section ^ how near the extremal horizon, the geometry is AdS2 x 
M 3 x T ' up to slowly varying logarithmic factors. We can thus ask what the effective 
dimensions of the operators dual to the modes given in table || are. Changing variables to 
the AdS2 coordinate y defined in (|5.9|), we see that 5X behaves for the six solutions as 



y 2 , y°, y, y 2 , and y + 2 , respectively. 
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Solution 


a 


bi 


b 2 


h 


Am 


I 

II 
III 

IV 
V 
VI 


5 

36 
5 

72 

5 

72 


5(HV5) 

144 
5(1-^5) 

144 


83 

30 
163 
60 

21 

20 



63-7^ 

120 
63+7^ 

120 


h + i 
h + i 
h 

bi 
bi 
h 


bx + l 
&i + l 

bi 

bi 

h 

bi 


2, source 
2, VEV 
0, VEV 
0, source 
, source 
l+ / E , VEV 



Table 2: The coefficients of the perturbative expansion (7.13) and the IR dimensions of the corre- 
sponding operators. The solution IV is in fact an exact pure gauge mode for which $^ = = 
const. 



Using (|5.16| ), the dimensions Air corresponding to various perturbations are given 
in the last column of table Since solution IV is exact, has A = 0, and is pure gauge, 
we suspect that the Air = modes we are seeing correspond to the conserved charge 
operator in the effective quantum mechanics. This is consistent with the fact that mode 
III, which produces a VEV of this operator, is seen to correspond to non-vanishing Qr. 
In this paper we only study solutions with vanishing R-charge, so the charged modes with 
Air = are not allowed. The remaining dimensions we find, 2 and 1+ 2 V ^ , correspond to 
irrelevant operators from the point of view of the IR near-^4dS*2 theory. The sources for 
such operators correspond to modes that fall off near the horizon as y~ l or y 2 . Since 
the operators are irrelevant, we expect that inducing them in the IR theory will not destroy 
the near-conformal IR solution we find. 



8. Discussion 

In this paper we initiated studies of black hole solutions charged under baryonic symmetries. 
Such solutions are asymptotic to AdS x Y, and the baryonic U(1)b symmetries appear due 
to the non-trivial topology of the Einstein space Y. We discussed the type IIB example 
AdS^ x T ' in some detail, and have also set up a consistent truncation for M-theory 
on AdSi x Q ' > . Perhaps our most surprising finding is that the type IIB charged 3- 
brane solution develops, in the zero-temperature limit, a novel kind of near-horizon region, 
which is a warped product AdSi x M 3 x T 1 ' with warp factors that are logarithmic in 
the AdS radius. This supergravity solution is smooth because the logarithms decrease the 
curvature of the solution; in fact, all curvatures approach zero at the horizon. In this sense 
this solution is reminiscent of the UV region of another solution based on the conifold, with 
a topologically non-trivial 3- form flux turned on |22|, ^5fl . That warped deformed conifold 
solution was super symmetric and automatically stable. In the present case, where the only 
non- vanishing supergravity fields are the metric and the self-dual 5-form, the solution does 
not seem to preserve any supersymmetry, and its stability is a serious issue. We carried 
out some highly non-trivial stability checks for our solution. 

We have shown that the simplest objects charged under the baryonic U(1)b, namely 
the wrapped D3-branes, do not condense. This still leaves the possibility that one of the 
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neutral fields might cause an instability. Our solution preserves a certain Z2 symmetry 
and we have checked stability with respect to one of the modes odd under the Z2. This 
well-known mode, dual to the operator Tr(Aj^4j — BjBj), turns on the difference of the sizes 
of the two 2-spheres that is present in a small resolution of the conifold, and also mixes 
with the U(1)r gauge field. We leave further studies of stability for future work. We also 
note that we have encountered difficulties in extending our numerical solution all the way 
to zero temperature. At the lowest temperature we have been able to reach numerically, 
Wh/2 ~ 0.85, which means that the neai-AdS2 throat is only beginning to develop. It 
would be interesting to construct the full numerical T = solution that matches onto 
the near-horizon form that we found analytically. One of our goals is to study the fermion 
modes in this background, following the interesting results of [0, ||, |9j|lC||. There are various 
extensions of this work related to finding other black brane solutions with baryonic charges. 
We hope to present these results in the future. 



Note Added 



After this paper was completed, we learned of the interesting work [33 where another 
potential instability of charged brane backgrounds was suggested. Such an instability, 
called the "Fermi seasickness" in [j33|, is caused by nucleation of a space-time filling D- 
brane towards the AdS boundary (for earlier discussions of similar D-brane instabilities see 
12, p4|). 4 In the dual gauge theory, this corresponds to an instability with respect to the 



Coulomb branch, where certain mesonic operators develop vacuum expectation values. In 
the background we have studied, this instability can be seen from computing the potential 
for probe D3-branes filling the (t, x) directions. Numerical computations show that for 
temperatures greater than about 0.2// the charged black branes constructed in this paper 
are stable with respect to such D3-brane nucleation, while for lower temperatures they 
become metastable. For any nonzero temperature the D3-brane is attracted near the 
horizon, which means that there exists a potential barrier that for large N prevents brane 
nucleation to AdS infinity. 
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A. Small Charge Limit 

Approximate solutions can be found in the small Q limit. Defining the dimensionless 

4 We thank Eva Silverstein for suggesting to us that our background may suffer from this instability. 
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(A.l) 



parameter Q = QL/r^, the small Q expansion takes the form 
w = Q 2 5w^ + QW 4) + • • • , 

5 = ^(i-J|)+qV 2) + qV 4) + --- , 

*=§(4-^)+^ (3)+ --- . 

rj = Q 2 5r]W + Q 4 <5r/ 4) + • • • , 

x = Q 2 Sx i2) + <?V 4) + • • • , 

where the starting point of the expansion obtained by setting Q = corresponds to the 
AdS-Schwarzschild solution. The small Q approximation (A.l) can be alternatively thought 
of as a large temperature expansion. 

To second order in Q, defining p = r/r^, one finds the solution 



Ju> (2) = , 5g {2) 



12r 4 



(i-p 2 ) . 



(2) _ _ 1 fl + 3p 2 -6p 4 



(A.2) 



40 



+ 3(1 -2/) log 



l + p 2 



This solution obeys the boundary conditions ( |4. 2[) at the conformal boundary and is also 
regular at the black hole horizon. 

The thermodynamic quantities fl4.4| ) become in this limit 



P 



Q 



2k\L 2 



Q 



4u> 



12 + Q< 



2AnL 2 



(24 - Q ; 



(A.3) 



One can check explicitly that the relation ([O]) is satisfied. It is also useful to note that 
the values rjh and Xh of rj and x at the horizon are 



Vh 



"" 3 q 2 



Xh 



10,8-2^ 



40 " /v 40 

B. A Consistent Truncation for Baryonic 2-Branes in AdS^ x Q 1 ' 1 ' 1 



(A.4) 



The consistent truncation of lid SUGRA that contains one of the two U(1)b gauge fields 
can be obtained as follows. The eleven-dimensional metric is 



ds 2 = e- 7 x/ 2 ds 2 M + 4L 2 e* 

pV -2X 



+ - 



8 

[del + sin 2 6zdcj>l) + 



\ + sin 2 6>i# 2 + del + sin 2 6> 2 # 2 ) 

x 2, (B.1) 



-675 



16 



COS I 



i=l 
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Here, ds 2 M is the line element on a non-compact four-dimensional space M. When x = V 



ll.l.l 



A = 0, (gj) describes M x Q 1 ' 1 ' 1 ||. 

The four-form F4 can be constructed using the forms W2 and ^5 defined as 

u>2 = - (sin#ic£#i A d^i — sin 6*2^2 A d^) , 

2 (B.2) 

o;5 = (dip + cos 9\d(j)\ + cos #2^2) A CJ2 A (sin ^3^3 A dcfo) . 
The four-form F4 is then 

F 4 = |e-f *vol M + ^e 2 ^-!^**/^) A w 2 . (B.3) 
Its dual, F7 = *F4, is given by 

3L 6 L 5 
F 7 = — uj 2 AW5 + A ^5 • (B.4) 

This ansatz preserves the Z2 space-time inversion symmetry where x^ — > —x^ is ac- 
companied by the interchange of two 2-spheres, (9±, <f>x) (62, fa)- This symmetry requires 
the metrics of these 2-spheres to be the same, thus eliminating an additional function. 



Equations (B.l) and ( |B,3| ) give a consistent truncation of eleven-dimensional SUGRA 
only when F A F = 0. The effective four-dimensional lagrangian is 

C eS = R- ie^-fri*, - |(^x) 2 - l^r?) 2 - 3(9,A) 2 - F(r?, X , A) , (B.5) 



where 



-8r,-2A-f x 



y(ry, X ,A) = 6 - " 2 9e 8 " +2A + 2e^ + e 6 ( A +^-16e 7 " +A+ ^-8e 7 ''+ 4A+ ^ . (B.6) 

One can again look for extrema of flB.5| ) of the form ( |3.1[ ), except now there are only 
two coordinates along the brane, dx 1 and dx 2 . The equations of motion resulting from C e g 
are 

„ ,/2 c,' 21r /2 63r /2 3r /2 \ $' 2 _ 2r) 3 V+2A 1 
x » + X '( 2 - + 1 + HIV + ^X' 2 + M - - 



r j 4 ' 16 A 2 y 21j 63^ 

x« + v (H + 1 + «v + -x' 2 + -a' 2 ) + ^W§*+« _ !«: 

\r g 4 ' 16 2 / 65 65 <9A 

. /I 21r /2 63r /2 3r /2 \ r<3?' 2 , 2 „_ 3 V , 2A r 
» + 9[~ + —V' 2 + ~^X 2 + yA' 2 J + — e -+ 2 " 2 *+ 2A + _y 

*" + *' (~ + I™' + 2V - |x' + 2A' 
Vr 2 '2 



" , + 2r (T r? ' 2 + il x/2 + ^ /2 



, 


(B.7a) 


, 


(B.7b) 


, 


(B.7c) 


, 


(B.7d) 


, 


(B.7e) 


. 


(B.7f) 
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Again, we can immediately integrate the equation, yielding 



(B.7g) 



where Q is an integration constant related to the charge of the black hole. Plugging ( |B.7g| ) 
into equations ( |B.7a|) - (|B.7d| ), we again eliminate all dependence on w, 



\r g 4 16 2 y Zlr^g 21g or) 



\r g 4 16 2 



Q 2 



-2r?+fx-2A 



A" + A' ( - + - + _ ?? - + 

r g i 15 



21r 4 # 



i : d' , 21r /2 , 63r /2 , 3r x/2 ^ ( Q 2 ^ 2 ??+|x-2A 



6r 4 <7 



4 

63c/ 9x 
\_dV_ 
6g d\ 



, 
, 
, 



(B. 



, / 1 21r /9 63r , 9 3r w9 \ Q 2 3 V 9> r 
5 +5 ( - + — V' 2 + ^-X' 2 + yA' 2 J + 2 " +5X 2A + ~^ = . 



r 4 ' 16 2 y 4H 2 

We hope to report on solutions of these equations in a later publication. 
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